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Abstract 

The Laplacian in an unbounded tubular neighbourhood of a hyperplane with non- 
Hermitian complex-symmetric Robin-type boundary conditions is investigated in 
the limit when the width of the neighbourhood diminishes. We show that the 
(~ l ' Laplacian converges in a norm resolvent sense to a self-adjoint Schrodinger opera- 

tor in the hyperplane whose potential is expressed solely in terms of the boundary 
coupling function. As a consequence, we are able to explain some peculiar spec- 
tral properties of the non-Hermitian Laplacian by known results for Schrodinger 
operators. 



1 Introduction 

, There has been a growing interest in spectral properties of differential operators in 

' shrinking tubular neighbourhoods of submanifolds of Riemannian manifolds, subject 

, to various boundary conditions. This is partly motivated by the enormous progress in 

' semiconductor physics, where it is reasonable to try to model a complicated quantum 

Hamiltonian in a thin nanostructure by an effective operator in a lower dimensional 
substrate. But the problem is interesting from the purely mathematical point of view 
as well, because one deals with a singular limit and it is not always obvious how the in- 
formation about the geometry and boundary conditions arc transformed into coefficients 
}J] , of the effective Hamiltonian. 

The interest has mainly focused on self-adjoint problems, namely on the Laplacian in 
the tubular neighbourhoods with uniform boundary conditions of Dirichlet [6l [5l 13 [17] 
or Neumann [H] or a combination of those [TS]. For more references see the review 
article The purpose of the present paper is to show that one may obtain an 

interesting self-adjoint effective operator in the singular limit even if the initial operator 
is not Hermitian and the geometry is rather trivial. 

We consider an operator iJ^ which acts as the Laplacian in a d-dimensional layer: 

H,u = -Au in M'*"^ x (0,e) , (1.1) 

where d ^ 2 and e is a small positive parameter, subjected to non-Hcrmitian boundary 
conditions on dfl^. Instead of considering the general problem, we rather restrict to a 
special case of separated Robin-type boundary conditions 

du 

— +ia{x')u = on dQ^ , (1-2) 



1 



where x — {xi, . . . ,Xd^i,Xd) = {x',Xd) denotes a generic point in and a is a real- 
valued bounded function. More precisely, we consider as the m-sectorial operator TJ^ 
on L2{^e) which acts as (jl.ip in the distributional sense on the domain consisting of 
functions u from the Sobolev space (fie) satisfying the boundary conditions (|1.2p . 
We postpone the formal definition to the following section. 

The model iJg in d = 2 was introduced in [3] by the present authors. In that paper, 
we developed a perturbation theory to study spectral properties of Hg with e fixed in 
the regime when a represents a small and local perturbation of constant ao (see below 
for the discussion of some of the results). Additional spectral properties of were 
further studied in |16j by numerical methods. The present paper can be viewed as 
an addendum by keeping a (and d) arbitrary but sending rather the layer width s to 
zero. We believe that the present convergence results provide a valuable insight into 
the spectral phenomena observed in the two previous papers. 

The particular feature of the choice (jl.2p is that the boundary conditions are 77- 
symmetric in the sense that H;, commutes with the product operator TT. Here CP denotes 
the parity (space) reversal operator (Tu){x) := u{x', £—Xd) and 7 stands for the complex 
conjugation (Tm)(x) := u{x); the latter can be understood as the time reversal operator 
in the framework of quantum mechanics. The relevance of non-Hermitian TT-symmetric 
models in physics has been discussed in many papers recently, see the review articles 
m [18] . Non-Hermitian boundary conditions of the type (|1.2p were considered in [13] 
to model open (dissipative) quantum systems. The role of (jl.2[) with constant a in the 
context of perfect-transmission scattering in quantum mechanics is discussed in |12| . 

Another feature of (|1.2p is that the spectrum of "does not explode" as the layer 
shrinks, meaning precisely that the resolvent operator (77^ -I- 1)~^ admits a non-trivial 
limit in L2{^e) as e — >■ 0. As a matter of fact, it is the objective of the present paper to 
show that converges in a norm resolvent sense to the (d — l)-dimensional operator 

Ho:=-A + a^ on L2iW^^^) , (1.3) 

which is a self-adjoint operator on the domain Wf (M'^"^). Again, we postpone the pre- 
cise statement of the convergence, which has to take into account that the operators 
and Hq act on different Hilbert spaces, till the following section (c/ Theorem 12. ip . 
However, let us comment on spectral consequences of the result already now. 

First of all, we observe that a significantly non-self-adjoint operator converges, in 
the norm resolvent sense, to a self-adjoint Schrodinger operator Hq. The latter contains 
the information about the non-self-adjoint boundary conditions of the former in a simple 
potential term. It follows from general facts [14] Sec. IV. 3. 5] that discrete eigenvalues 
of Hg either converge to discrete eigenvalues of Hq or go to complex infinity or to the 
essential spectrum oi Hq as e — >■ oo. 

In particular, assuming that i/^ and Hq have the same essential spectrum (indepen- 
dent as a set of e), the spectrum of H^ must approach the real axis (or go to complex 
infinity) in the limit as e — )■ 0. Although numerical computations performed in |16| sug- 
gest that Hg might have complex spectra in general, perturbation analysis developed 
in [3] for the 2-dimensional case shows that both the essential spectrum and weakly 
coupled eigenvalues are real. The present paper demonstrates that the spectrum is real 
also as the layer becomes infinitesimally thin, for every d ^ 2. We would like to stress 
that the TT-symmetry itself is not sufficient to ensure the reality of the spectrum and 
that the proof that a non-self-adjoint operator has a real spectrum is a difficult task. 

It is also worth noticing that the limiting operator Hq provides quite precise infor- 
mation about the spectrum of Hg in the weak coupling regime for d — 2 and \a\ < ir/e 
(e fixed). Indeed, let us consider the following special profile of the boundary function: 

adx') := Qo + c(3{x') , 

where ao is a real constant, ^ is a real- valued function of compact support and c is a real 
parameter (the regime of weak coupling corresponds to small c) . Note that the essential 
spectrum of both Hg and Hq coincides with the interval [q;§,oo), for /3 is compactly 
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supported (c/ [H Thm. 2.2]). It is proved in g] that if ao J^l3{x')dx' is negative, 
then possesses exactly one discrete real eigenvalue ^(c) converging to a§ as c — >■ 0+ 
and the asymptotic expansion 

fi{c) ^ al~ c^al(^J P{x')dx'^ + O(c^) as c 0+ 

holds true. As a converse result, it is proved in g] that there is no such a weakly coupled 
eigenvalue if the quantity ag Jj^l3(x') dx' is positive. These weak coupling properties, 
including the asymptotics above, are well known for the Schrodinger operator Hq with 
the potential given by a^, see [5]. 

At the same time, the form of the potential in Hq explains some of the peculiar 
characteristics of iJ^ even for large c. As an example, let us recall that a highly non- 
monotone dependence of the eigenvalues of on the coupling parameter c was observed 
in the numerical analysis of |16| . As the parameter increases, a real eigenvalue typically 
emerges from the essential spectrum, reaches a minimum and then comes back to the 
essential spectrum again. This behaviour is now easy to understand from the non-linear 
dependence of the potential on c. 

On the other hand, we cannot expect that Hq represents a good approximation of H^ 
for the values of parameters for which H^ is known to possess complex eigenvalues |16j . 
It would be then desirable to compute the next to leading term in the asymptotic 
expansion oi H^ as e — >■ 0. 

This paper is organized as follows. In Section [2] wc give a precise definition of the 
operators TJg and Hq and state the norm resolvent convergence of the former to the 
latter as e fThcorcm l2.ip . The rest of the paper consists of Section |2] in which a 
proof of the convergence result is given. 

2 The main result 

We start with giving a precise definition of the operators H^ and Hq. 

The limiting operator (|1.3p can be immediately introduced as a bounded perturba- 
tion of the free Hamiltonian on L2{R'^^^), which is well known to be self-adjoint on the 
domain W^iR''-^). For later purposes, however, we equivalently understand Hq as the 
operator associated on L2(M'*~^) with the quadratic form 

hQ[v]:= f \V'v{x')\^dx' + ( a{x'f\v{x')\dx', 

V e £>(/io) := W^iR'^-^). 

Here and in the sequel we denote by V' the gradient operator in M''^^, while V stands 
for the "full" gradient in R'^. 

In the same manner, we introduce i/g as the m-sectorial operator associated on 
i2(^^e) with the quadratic form 

h^[u] / \\/ u{x)\'^ dx + i / a{x)\u{x\e)\'^dx'~i / a{x )\u{x',0)\'^dx', 
ue^{hQ) :^W^{n,). 

Here the boundary terms are understood in the sense of traces. Note that TJg is not 
self-adjoint unless a = (in this case iJg coincides with the Neumann Laplacian in 
the layer fie). The adjoint of H^ is determined by simply changing a to —a (or i to 
— i) in the definition of h^. Moreover, if^ is T-self-adjoint [TJ Sec. HI. 5] (or complex- 
symmetric [lO]), i.e. H* = TiJeT. 

The form he is well defined under the mere condition that a is bounded. However, if 
we strengthen the regularity to a € VF^(R'^~^), it can be shown by standard procedures 
(c/ [H Sec. 3]) that H^ coincides with the operator described in the introduction, i.e.. 
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it acts as the (distributional) Laplacian on the domain formed by the functions u 
from VF|(Oe) satisfying the boundary conditions (|1.2p in the sense of traces. 

The operator Hq is clearly non-negative. An analogous property for is contained 
in the following result 

cr(iJ£) C |z e C I Rez ^ 0, |Imz| < 2|la|loo%/ReT| . (2.1) 

Here and in the sequel we denote by || • ||oo the supremum norm. (j2.1[) can be proved 
exactly in the same way as in [?J Lem. 3.1] for d = 2 by estimating the numerical range 
of H^. In particular, the open left half-plane of C belongs to the resolvent set of both 
and iJo- 

Another general spectral property of H^, common with iJo, is that its residual 
spectrum is empty. This is a consequence of the T-self-adjointness property of as 
pointed out in [U Corol. 2.1]. 

Since and Hq act on different Hilbert spaces, we need to explain how the con- 
vergence of the corresponding resolvent operators is understood. We decompose our 
Hilbert space into an orthogonal sum 

L2{^e)=S--)e®S)^ , (2.2) 

where the subspace S)^ consists of functions from L2{^e) of the form x i— >■ 'ip{x'), i.e. 
independent of the "transverse" variable Xd- The corresponding projection is given by 

€ 

{P,u){x) ~ ^ J u{x) dxd (2.3) 



and it can be viewed as a projection onto a constant function in the transverse variable. 
We also write Pj^ := I — P^. Since the functions from i^^ depend on the "longitudinal" 
variables x' only, i^^ can be naturally identified with L2(M'*^^). Hence, with an abuse 
of notations, we may identify any operator on L2(R^^^) as the operator acting on 
S)s C L2(f^e), and vice versa. 

The norm and the inner product in L2{^le) will be denoted by || • \\^ and {■,-)e, 
respectively. We keep the same notation || • \\^ for the operator norm on L2{^e)- The 
norm and the inner product in L2{M.'^~^) will be denoted by || • || and (•, •), respectively, 
i.e. without the subscript e. All the inner products are assumed to be linear in the first 
component. Finally, we denote the norm in W2{^e) by || • and we keep the same 
notation for the norm of bounded operators from L2{^e) to W2{^e)- 

Now we are in a position to formulate the main result of this paper. 

Theorem 2.1. Assume a G W^{M.'^^^) . Then the inequalities 

\\{H, + l)-i - {Ho + l)-'Pe\\^ ^ Ce, (2.4) 
||(i/, + l)-i-(l + Q)(Fo + l)-'PeL_i ^C(£)e (2.5) 

hold true, where Q{x) := -~ia{x') Xd and 



l_ , (||V'a||oo + 2||a||oo)2 

r2 




Ci(e) 
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Let us discuss the result of this theorem. It says that the operator 77^ converges 
to Hq in the norm resolvent sense. Note that, contrary to what happens for instance 
in the case of uniform Dirichlet boundary conditions, here we can choose the spectral 
parameter fixed {e.g. —1 S p{H^) D p{Ho) as in the theorem) and still get a non-trivial 
result. 

If we treat the convergence of the resolvents in the topology of bounded operators 
in L2{^e), the estimate (|2.4p says that the rate of the convergence is of order 0(e). At 
the same time, if we consider the convergence as for the operators acting from L2{^e) 
into W2{fle), to keep the same rate of the convergence, one has to use the function Q. 
This functions is to be understood as a corrector needed to have the convergence in 
a stronger norm. Such situation is well-known and it often happens for singularly 
perturbed problems, especially in the homogenization theory, see, e.g., [51 [^. 

3 Proof of Theorem D 

Throughout this section we assume a € Wl^{R'^~^). With an abuse of notation, we 
denote by the same symbol a both the function on M'*^^ and its natural extension 
X ^ a{x') to WK 

We start with two auxiliary lemmata. The first tells us that the subspace i^^ is 
negligible for in the limit as e — > 0. 

Lemma 3.1. For any f e L2{^e), we have 

\\{H, + ir'ptf\i,^-\\ptih. (3.1) 

Proof. For any fixed / e L2(f^e), let us set u := {H, + ly^P^J £ S)(i?e) C W^^in^). 
In other words, u satisfies the resolvent equation 

V« e Wl{n,) , h,{u,v) + (u,«), = {Ptf.v), , 

where /ie(-,-) denotes the sesquilinear form associated with the quadratic form /ie[-]. 
Choosing u for the test function v and taking the real part of the obtained identity, we 
get 

\\u\\i^ = Jic{P^f,u)^^Rc{P^f,P^u)^ ^ \\P^fUP^u\\,. (3.2) 

Employing the decomposition u = P^u + P^u, the left hand side of (|3.2p can be esti- 
mated as follows 

Ml, > W^Al ^ \\d,ut = WddPM'e > {^lef \\PM\1 ■ (3.3) 

Here the last inequality follows from the variational characterization of the second eigen- 
value of the Neumann Laplacian on L2{{Q,e)) and Fubini's theorem. Combining (|3.3p 
with (|3.2p . we obtain 

\\P,UU^{e/7rr\\P,^f\U. 

Finally, applying the obtained inequality to the right hand side of p.2[) . we conclude 
with 

\\u\\l,^ie/nr\\P,^f\\l. 
This is equivalent to the estimate p.ip . □ 

In the second lemma we collect some elementary estimates we shall need later on. 
Lemma 3.2. We have 

|e---^ -IK (3.4) 
\e-''^^--l + iax,\^^\\a\\lxl, (3.5) 
|V(e-'"^'' -l + iaxrf)| ^ ||a||ooa;dy^!|Q||^ + ||V'a||^a;2. (3.6) 
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Proof. The estimates p.4p and (j33j) arc elementary and we leave the proofs to the 
reader. The last estimate (|3.6p follows from p.4p and the identity 



y^g-iax^ _ ^ _^ ^^^^^ ^ _ g-io^d) 



taken into account. □ 

Wc continue with the proof of Theorem 12.11 Let / e L2{^e). Accordingly to (|2.2|) . 
/ admits the decomposition f = Pef + f and we have 

11/11? = WPJWl + WP^fWl = eWP^ff + WPtfWl . (3.7) 

We define u := {Hi. + 1)^^/ and make the decomposition 

u^uo + ui with uo := {H, + 1)-^ PJ , m := [H, + ly^ P^ f . (3.8) 

In view of Lemma l3. 11 ui is negligible in the limit as e — > 0, 

\\u4,.,^-\\P,^f\U. (3.9) 

TT 

It remains to study the dependence of uq on e. We construct uq as follows 

=e-'"(^')^''wo(x') + 1^1(3;), where wq := {Hq + ly^ P,f (3.10) 

and wi is a function defined by this decomposition. 

First, we establish a rather elementary bound for wq. 

Lemma 3.3. We have 

Ikolle,! ^ WPelWe- 

Proof. By definition, ifo satisfies the resolvent equation 

■iv^.W^iM."-^), hoiwo,v) + iwo,v)^iPJ,v), (3.11) 

where ho{-,-) denotes the sesquilinear form associated with the quadratic form ho[-]. 
Choosing wq for the test function w, we get 

\\Vwof + \\awof + \\wof = (Pef, wo) s; \\Pef\\\\wo\\ ■ (3.12) 

In particular, 

W^oW ^ WPefW- 

Using this estimate in the right hand side of p.l2p . we get 

\\V'wof + \\wof^\\PJ\\\ 
Reintegrating this inequality over (0, e), we conclude with the desired bound in il^. □ 

It is more difficult to get a bound for wi . 
Lemma 3.4. We have 



\wi\yi^C\e\\PJ\U with Co:= 



IV'alloo 



V3 

Proof. By definition, uq satisfies the resolvent equation 

yveW^in,), h,iuo,v) + {uo,v)e = [Pef.v),. 

Choosing wi for the test function v and using the decomposition p.lOp . we get 

he[wi\ + \\wi\\l = {Pef,Wi)e - K{uo - Wi,Wi) - {uq - Wi,Wi)e =■ . (3.13) 



It is straightforward to check that 
with 

F^ ■.= i{xdWoVa,e"'^^Vwi)^-i{Vwo,Xde'''^'wiVa)^. 

Here the first equahty follows by algebraic manipulations using an integration by parts, 
while the second is a consequence of p. lip , with x' i— >■ e'"*^^ '^^'''Wi{x' ,Xd) being the test 
function, and Fubini's theorem. At the same time, {uq — wi,wi)s = (wq, e'"^''wi)£. 
Hence, 

F,^F^+{PJ,w,-e"'^''wi)^. 
We proceed with estimating F^: 

3 /2 

^ ^(|lV'a|loc|ko!ll|V'i.i|U + \\V'a\U\V'wo\\\\w^\U + ||a|U|lP,/|l||^«i|l, 
< ^(||V'a|UVll«^olP + ||V'«;o|| Vl|V'«^i||2 + \\w,\\^, + \\a\U\Pef\\h'i\\e) 



^ -^(||V'a||oo||^«o||e,l + ll"lloo||Pe/||e)||u'l||e,l. 

Here the first inequality follows by the Schwarz inequality, an explicit value of the 
integral of x'^ and obvious bounds such as p.4p . 

Finally, taking the real part of (|3.13p and using the above estimate of \F^\, we get 



1 S; ^(\\Va\\^\\wo\\eS + \Hoc.\\Pef\\ 



The desired bound then follows by estimating HuioUe,! by means of Lemma 13.31 □ 

Now we arc in a position to conclude the proof of Theorem l2.1l bv simply comparing u 
with wq. As for the convergence in the topology of L2{^le), we write 

\\U - WoWe = \\ui+Wi + {C'"'-' ~ iKIle 

^ \\u,\U + \\wi\\, + \\{e-'"''^-l)wo\\e. 



Here the last term can be estimated using Lemma [3^21 as follows 

3 /2 

||(e---^ - iVolle ^ ^ llalUlkoll = ^ llalloolkol 

Hence, using p.9p . Lemma [33] and Lemma [3. 31 we get the bound 

s 

\\u -WoWe ^ ll^tllle.l + II Willed + --j= ||a|| oo || || e 



^ e lln^/ll. + ^(llV'alloo + 2||a||oo) ||P./||e) 
^Ce\\f\U 



Here the last estimate follows by the Schwarz inequality recalling l\'3.7\i and holds with 
the constant C as defined in Theorem 12. II This proves (|2.4p . 
As for the bound (12.51). wc have 



||it - (1 + Q)wo\\eA = \\ui +wi + {e - 1 + iaXd)Wa\\e,l 

||ui||e,i + ||wi||c,i + IKe"'"^" - 1 + iaxd)wo\\,^i . 
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Here the last term can be estimated using Lemma 13.21 as follows. Employing the indi- 
vidual estimates 

\\{e-'^^^ - 1 + iax,)wo\U ^ 0^ llallLll'^oll = ^ ll«llLll^o||s , 

11 (e--<^ - 1 + ,ax,)Vwo\U ^ ^ V'u;o|| = ^ V'u;o||e , 

=■3/2 



|l«;oV(e---'' - 1 + iax,)\U ^ \\a\\^VML+WW^ \\wo\\ 



= ^||a||^Vll«ll^ + l|V'a||^£2||«;o||e, 
and the Schwarz inequality, we may write 

- 1 + iaxd)wo\U.i ^ Ci(e) e \\wo\U,i 

with the same constant Ci (e) as defined in Theorem 12.11 Consequently, using 
Lemma 13.41 Lemma |3.3l and the Schwarz inequality employing (j3.7p . we get the bound 

\\u-{l + Q)wo\U,i^C{e)e\\f\U 

with 

Cie) y^ + (Co + Ci(e))' . (3.14) 

Note that C(e) coincides with the corresponding constant of Theorem 12.11 This con- 
cludes the proof of Theorem 12. II 
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